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Abstract 

In this work we present a formal generalization of the Hamilton- Jacobi formalism, 
recently developed for singular systems, to include the case of Lagrangians con- 
taining variables which are elements of Berezin algebra. We derive the Hamilton- 
Jacobi equation for such systems, analizing the singular case in order to obtain 
the equations of motion as total differential equations and study the integrability 
conditions for such equations. An example is solved using both Hamilton-Jacobi 
and Dirac's Hamiltonian formalisms and the results are compared. 



1 Introduction 

In this work we intend to study singular systems with Lagrangians containing 
elements of Berezin algebra from the point of view of the Hamilton- Jacobi for- 
malism recently developed [|I], |^. The study of such systems through Dirac's 
generalized Hamiltonian formalism has already been extensively developed in lit- 
erature [^, ^, ^ and will be used for comparative purposes. 

Despite the success of Dirac's approach in studying singular systems, which 
is demonstrated by the wide number of physical systems to which this formal- 
ism has been applied, it is always instructive to study singular systems through 
other formalisms, since different procedures will provide different views for the 
same problems, even for nonsingular systems. The Hamilton- Jacobi formalism 
that we study in this work has been already generalized to higher order singular 
systems [||, |^ and applied only to a few number of physical examples as the elec- 
tromagnetic field 1^, relativistic particle in an external Electromagnetic field 
and Podolsky's Electrodynamics ^. But a better understanding of this approach 
utility in the studying singular systems is still lacking, and such understanding 
can only be achieved through its application to other interesting physical systems. 

Besides that, Berezin algebra is a useful way to deal simultaneously with 
bosonic and fermionic variables in a unique and compact notation, what justifies 
the interest in studying systems composed by its elements using new formalisms. 

The aim of this work is not only to generalize the Hamilton- Jacobi approach 
for singular systems to the case of Lagrangians containing Berezinian variables 
but to present an example of its application to a well known physical system, 
comparing the results to those obtained through Dirac's method. 



We will start in section 2 with some basic definitions and next, in section 
3, we will introduce the Hamilton- Jacobi formalism to Berezinian systems using 
Caratheodory's equivalent Lagrangians method. In section 4 the singular case 
is considered and the equations of motion are obtained as a system of total dif- 
ferential equations whose integrability conditions are analyzed in section 5. The 
equivalence among these integrability conditions and Dirac's consistency condi- 
tions will be discussed separately in the appendix. In section 6 we present, as 
an example, the electromagnetic field coupled to a spinor, which is studied using 
both the formalism presented in this work and Dirac's Hamiltonian one. Finally, 
the conclusions are presented in section 7. 

2 Basic definitions 

We will start from a Lagrangian L (q, Qj that must be an even function of a set 
of A^ variables g* that are elements of Berezin algebra. For a basic introduction 
in such algebra we suggest the reader to refer to ref. 0], appendix D, from which 
we took the definitions used in this paper. A more complete treatment can be 
found in ref. [ITU . 



The Lagrangian equations of motion can be obtained through variational prin- 
ciples from the action S = J Ldt 

6rS drL d drL 



6q' dq^ dt Q q 

were we must call attention to the use of right derivatives. 

The passage to Hamiltonian formalism is made, as usual, by defining the 
momenta variables through right derivatives as 

Pi = —i (2) 

dq 



and introducing the Hamiltonian function as (summing over repeated indexes) 

H = p,q -L, (3) 

were the ordering of momenta to the left of velocities shall be observed since they 
were defined as right derivatives. This ordering will be, of course, irrelevant when 
we deal with even elements of the Berezin algebra. The Hamiltonian equations 
of motion will be given by 

■ i diH -.-.p,, drH . drH -.-.p,, diH 

If we use the Poisson bracket in the Berezin algebra, given by 

we get the known expressions 

q=[q\H]^- P,= {p,,H}^. (6) 

For simplicity and clarity we will refer to these brackets as Berezin brackets. 
These brackets have similar properties to the usual Poisson brackets 

{F,G'}^ = -(-lf(-)^(«){G,F}^, (7) 

{F, GK}^ = {F, G}^ K + (-l)^(-)^(«) G {F, K}^ , (8) 

+ (-lf(--)^(«){K,{F,G'}^}^ = 0, (9) 

were the last expression is the analogue of Jacobi's identity. 

Similarly to the usual case, the transition to phase space is only possible if the 
momenta variables, given by Eq. (0), are independent variables among themselves 



so that we can express all velocities Qi as functions of canonical variables (?\pj). 
Such necessity implies that the Hessian supermatrix 

^. - % - -^ (10) 

dq' dq dq' 

must be nonsingular. Otherwise, if the Hessian has a rank P = N — R, there will 
be R relations among the momenta variables and coordinates g* that are primary 
constraints (that we suppose to have definite parity), while R velocities q will 
remain arbitrary variables in the theory. The development of Dirac's General- 
ized Hamiltonian Formalism is straightforward: the primary constraints have to 
be added to the Hamiltonian, we have to work out the consistency conditions, 
separate the constraints in first and second-class ones and define the Dirac brack- 
ets using the supermatrix whose elements are the Poisson brackets among the 
second-class constraints 0]. 

3 Hamilton-Jacobi formalism 



From Caratheodory's equivalent Lagrangians method [^ we can obtain the Ha- 
milton-Jacobi equation for the even Lagrangian L (q,qj. The procedure is similar 
to the one applied to usual variables: given a Lagrangian L, we can obtain a 
completely equivalent one given by 

were S (g*, t) is an even function in order to keep the equivalent Lagrangian even. 
These Lagrangians are equivalent because their respective action integrals 
have simultaneous extremum. Then we choose the function S (g*, t) in such a way 
that we get an extremum of L' and, consequently, an extremum of the Lagrangian 
L. 



For this purpose, it is enough to find a set of functions f3\q^ , t) and S {q\ t) 
such that 

L'{q\q = (3\q^,t),t)=0 (11) 

and for all neighborhood oi q = (3^{q-^ , t) 



L'{q\q)>0. 



(12) 



With these conditions satisfied, the Lagrangian L' (and consequently L) will 
have a minimum in ^ = /3'^{q\ t) so that the solutions of the differential equations 
given by 

will correspond to an extremum of the action integral. 
From the definition of L' we have 

dS{q^,t) drS{q^,t)dq' 



L' = L(q\q') 



dt dqi dt ' 

where again we must call attention to the use of the right derivative. 
Using condition ([TT]) we have 



L{q^,q') 



dS{q\t) drS{q^,t) . 



q 



dS_ 
'dt 



q =f3i 



dt dq' 

drS{q\t) 



0, 



q =f3i 



L{q',q') 



dq^ 



q =l3i 



In addition, since L' has a minimum in 9 = /?*, we must have 



drL' 



dq 



q =l3i 



drL dr ( dS^ 



.dq dq\dtj. 



0, 



q =/3' 



or 



drL drS{q^,t) 

Id q' dq' 

drS{q\t) 



dq^ 



q =13' 



q =p 

drL 
dq' 



q =13' 



(13) 



(14) 
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Now, using the definitions for tlie conjugated momenta given by Eq. @, we 






get 



We can see from tliis result and from Eq. ([T^) that, in order to obtaining an 
extremum of the action, we must get a function S {q^ , t) such that 

dS , , 

a = -^» <!«' 

where ffo is the canonical Hamiltonian 

Ho = Piq -L(q\q') (17) 



and the momenta Pi are given by Eq. (p^. Besides, Eq. (|T6D is the Hamilton- 
Jacobi partial differential equation (HJPDE). 

4 The singular case 

We now consider the case of a system with a singular Lagrangian. When the 
Hessian supermatrix is singular with a rank P = N—Rwe can define the variables 
g* in such order that the P x P supermatrix in the right bottom corner of the 
Hessian supermatrix be nonsingular, i.e. 

det \Hat\ + 0; E,^ = .; ., ; a, 6 = i? + 1, ..., N. (18) 

dq dq 

This allows us to solve the velocities q as functions of coordinates q's and 
momenta Pa, i.e., q = f"" {q\Pb)- 

There will remain R momenta variables Pa dependent upon the other canonical 
variables, and we can always 0, §, |12[ write expressions like 



Pa = -Ha (q'; Pa); " = 1, •••, R] (19) 



that correspond to the Dirac's primary constraints ^a = Pa + Ha (g*; Pa) ~ 0. 
The Hamiltonian Ho, given by Eq. (|17D, becomes 

Ho = Par + Pa\,^=-H, ■ '?" -L {q\ Qa, %= P) , (20) 

where a, /? = 1, ..., i?; a = R+ 1, ..., N. On the other hand we have 

drHo _ drt , _ dr£ _ drl^drf^ _ 

dq^ '^^dq'^^^" dq" d q" d q" ~ ' 
so the Hamiltonian Ho does not depend exphcitly upon the velocities cfa- 

Now we will adopt the following notation: the time parameter t will be called 
t^ = q^; the coordinates g" will be called t"; the momenta Pa will be called Pa 
and the momentum po = Pq will be defined as 

Then, to get an extremum of the action integral, we must find a function 
S [ta] q"', t) that satisfies the following set of HJPDE 

H',^Po + Ho U t^; q^; P- = ^= 0' (22) 

H'a^Pa + Ha U t^', 9^ P- = ^] = 0" (23) 

where a, jS = 1, ..., R. If we let the indexes a and /? run from to i? we can write 
both equations as 

H'a^Pa + Ha (t^; ^^ Pa = ^^ = 0- (24) 

From the above definition and Eq. (^0|) we have 

Oj-iif) Oy-LOj-f ^-.Pru-.P, \ (JrHa -a Or t , ^ < P,, , -h 

' 'J — •' _ (^_]^^M6)M«) _i q _|_p -L (—\\w q 



dpb d (f dpb dpb " dp, 



'b 



dpb dpb 
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where we came back to a = 1 , . . . , i?. 

Multiplying this equation by dt = dt^ and (—1) ''', we have 

OPb Opb 

Using t" = g°, letting the index a run again from to i? and considering 

P(o) = -P(tO) = we have 

dq' = (-l)^w+^('')^w ^de. (25) 

OPb 

Noticing that we have the expressions 

dqf' = (-l)^(«+^(«^«^) M^dr = (-l)^(«+^(«^('^) d^dr = dt'' 

opp 

identically satisfied for a, /3 = 0, 1, ..., i?, we can write Eq. (^) as 

dq' = (_l)^«+^w^w ^^de- z = 0, 1, ..., A^. (26) 

OPi 

If we consider that we have a solution S {q\t) of the HJPDE given by Eq. 



( p4D then, differentiating that equation with respect to g*, we obtain 



a?* aP^ dtf^dq' dpa dq^dq^ ^ ' 

for a, /3 = 0, 1, ..., R. 

From the momenta definitions we can obtain 

Now, contracting equation (^) with dt" (from the right), multiplying by 
{^—V) '*' '"' and adding the result to equation ( P8| ) we get 



dp, + (-1)^«^('^) ^-^dt" 
ag* 



O^b I , „ , _xP, .^P, , , , X ^P, ,^P, , Vp..,4^P, ,"1 / _nP, .^P. , Orl2„, 



dq" - (_i)^w-P(-) (_i)(^M+^w)(^w+^(«)) (_l)^W^(a) ' "^ df^ 
dq'^dq"- \ dp a , 

d^ / . _ ._.^p„P(„) /_iUP(.)+P(^))(pw+P(^)) / 1^Pw■P(« ^lEkdr 
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dq^dq 






dq^dtP 
where we used the fact that 






(29) 



dqWq^ 
and that we have the following parities 

\dq^dq^ J 



^rS = (_ 1)^(0^0) _^^ 



dpj 



dqdq^ ' 



P« + P0-), 



^H + ^0> 



If the total differential equation given by Eq. ( PUD applies, the above equation 
becomes 

(30) 



^P(,^P|r.^ ^rH^ 



dp, = - (-l)-«-(") ^de; ^ = 0, 1, ..., iV. 

Making Z = S [t"'] q"') and using the momenta definitions together with Eq. 
we have 

dZ = ^dt^ + ^dq\ 



dt^ dq'' 



dZ = -Hpdt^ +pa I (-l)^w+^w^w "^^de ] . 
With a little change of indexes we get 



dpa 



dZ = {-Ha + (-l)^W+^W^(«p,M^^ df. 
\ dpa J 



(31) 



This equation together with Eq. (|26|) and Eq. (|30|) are the total differential 
equations for the characteristics curves of the HJPDE given by Eq. (^4]) and, 
if they form a completely integrable set, their simultaneous solutions determine 
S {t"'; g°) uniquely from the initial conditions. 
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Besides that, Eq. ( p6|) and Eq. (^) are the equations of motion of the system 
written as total differential equations. It is important to observe that, in the 
nonsingular case, we have only H'q j^ and no others H'^; so that these equations 
of motion will reduce naturally to the usual expressions given by Eq. (^ . 

5 Integrability conditions 

The analysis of integrability conditions of the total differential equations (p6D, 
( PPD and (^) can be carried out using standard techniques. This have already 
been made p|, [13|, [l^ for systems with usual variables, and here we will present 



the analysis of the integrability conditions for Berezinian singular systems. 
To a given set of total differential equations 

dq' = A\ {tf^, q^) dr (32) 

(i, j = 0, 1, ..., A^ and a, (3 = 0, 1, ..., R < N) we may associate a set of partial 
differential equations [^ 

XJ = ^A\ = 0, (33) 



where X^ are linear operators. 



Given any twice differentiable solution of the set (|33|), it should also satisfy 
the equation 

[X^,Xf,]f = 0, (34) 

where 

[X„, Xp] f = (x^Xp - (-1)^('^)^(« XpX^ f (35) 



is the bracket among the operators X^,. This implies that we should have 

[X^,Xp]f = CyX,f. (36) 
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So, the commutation relations ( ^5]) will give the maximal number of linearly 
independent equations. Any commutator that results in a expression that can't 
be written as Eq. ( |36D must be written as a new operator X and be joined to 
the original set (P^D, having all commutators with the other X's calculated. The 
process is repeated until all operators X satisfy Eq. (^). 

If all operators X^ satisfy the commutation relations given by Eq. (|36| ) the 
system of partial differential equations (^) is said to be complete and the corre- 
sponding system of total differential equations (|32|) is integrable if, and only if, 
the system ( |33D is complete. 

Now, we consider the system of differential equations obtained in the previous 
section. First we shall observe that if the total differential equations ( ]26|) and (|30D 



are integrable the solutions of Eq. (|3TD can be obtained by a quadrature, so we 
only need to analyze the integrability conditions for the last ones, since the former 
will be integrable as a consequence. 

The operators X^ corresponding to the system of total differential equations 
formed by Eq. ( PB] ) and Eq. (|3D|) are given by. 



Xj(tf,q\p- 



dq^ dpi dpi dq^ 

(9g* dpi 



(_l)no+^w-P(/) (_i)^wn^) (_i)(n.)+no)(^{/)+no) 



drK dif 



dq^ dpi ' 

Xj(t'^,q'',p'') = M.^lEk _ (_i)n.)n/) M^M, 
" ^ ' ' ^ dq^ dpi dq^ dpi ' 

X^f{tP,q^,p^)={f,H'J^; (37) 

where i = 0,1, ..., N; a, /3 = 0, 1, ..., i?; a = R + 1, ..., N and we have used Eq. 
( ^2|) and Eq. (|33D together with the result 

^ = (-!)''« (-1)^«^(-' 1^. (38) 

dq' ^ ' ^ ' dq' ^ ' 
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It is important to notice that the Berezin bracket in Eq. (37) is defined in a 
2N + 2 dimensional phase space, since we are including g'' = t as a "coordinate". 
Now, the integrability condition will be 

= X^ {/, H'p]^ - (-1)^(")^(« Xp {/, K}^ = 0, 

[X., X,] / = {{/, H',}^ , K}^ - (-1)^(")^(« {{/, H'J^ , i/;,}^ = 0, (39) 
that will reduce to 



[XaiX/s] f 



(_l)(^a)-.^(.))^(.){,:,^,{;,^;^}J^ 

(_l)^w^(/3) (_i)(%)+^w)n«+^(/)^w {iJ^;,!//;^,/}^}^ = 0, 



[Xa, x,]f = - (-i)^(/.^(") (_i)^(.)^M {i/;, {/, H',}^}^ 

_ (_l)^(<.)na) (_l)^(/)^(«+^«.)^(/3)+^(/)^«.) ^H'p, {H'^, /}^}^ = 0, 



Pin Pi, 



[X^,X^]f = -(-l)-(/)-w 



(_l)^(.,^,.,{i^;,{;,^;^}^} 



B 



B 



+ (-if(/'^(«{i/;,,{i7;,/}^}J=o, 
[x«,x,] / = - (-i)^(/)^(") [- (-i)^(/)^(") {/, [h'.^k]^}^ 

[X«, Xp] f = {/, {h'^, H'^}^}^ = 0, (40) 

when using the Jacobi relations for Berezin brackets given by Eq. (|^) and the 
fact that 

P{{A,B}^) = PiA) + PiB). (41) 

So, the integrability condition will be 



{h'^,K}^ = 0; \/a,p. 



(42) 
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It is important to notice that the above condition can be shown to be equiva- 
lent to the consistency conditions in Dirac's Hamiltonian formahsni but, to keep 
the continuity of the presentation, we will postpone the demonstration of this 
equivalence to the appendix. 

Now, the total differential for any function F (t^, q"',p"'] can be written as 

d F 8 F 8 F 

dF = ^dq^ + ^dp^ + ^de, (43) 



\ dq" dpa dq'' dp'' (9t° / 

dF = {F,H'J^de- (44) 

where the Berezin bracket above is the one defined in the 2N + 2 phase space 
used in Eq. (pT]). Using this result we have 

dH'p = {H'p,H'J^^de (45) 

and, consequently, the integrability condition (^2]) reduces to 

dH'^ = 0, Va. (46) 

If the above conditions are not identically satisfied we will have one of two 
different cases. First, we may have a new H' = 0, which has to satisfy a condition 
dH' = 0, and must be used in all equations. Otherwise we will have relations 
among the differentials df^ which also must be used in the remaining equations 
of the formalism. 
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6 Example 

As an example we analyze the case of the electromagnetic field coupled to a 
spinor, whose Hamiltonian formalism was analyzed in references |^, Q. We will 
consider the Lagrangian density written as 

L = —F^uF'^" + z V' t'^ (^M + ie\) ip-m^j^, (47) 

where A^ are even variables while ip and ifj are odd ones. The electromagnetic 
tensor is defined as F^^ = d'^A" — d'^A'^ and we are adopting the Minkowski 
metric rj^^, = diag{+l, —1, —1, —1). 

6.1 Hamiltonian formalism 

Let's first review Dirac's Hamiltonian formalism. The momenta variables conju- 
gated, respectively, to A^, ip and -0, are 

P, = Tr-r^ = -Fo,, (48) 

o A 

p, = §:^ = zV^/, P. = — = 0, (49) 

dip ■" 5^ 

where we must call attention to the necessity of being careful with the spinor 
indexes. Considering, as usual, V' as a column vector and -ip as a row vector 
implies that p^p will be a row vector while p- will be a column vector. 
From the momenta expressions we have the primary constraints 

(j)i=Po^ 0, (50) 

(p2=Pi,-ii^ 7°, 03 = P^ ~ 0. (51) 

The canonical Hamiltonian is given by 

He = Jncd'x = J (p, a' + (p^)^ W^ + ( P^)^ (^)^ - ^) d'x, (52) 
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He = jHcd'x = J^^F,,F'^-^p,p'-Ao[d,p'-ei:-f'i;) 

—i ^ 7"' (dj + ieAj) tp + m 'tp ijj d^x. (53) 

The primary Hamiltonian is 

Hp = J (He + Ai0i + </>2A2 + A303) d'x, (54) 

where Ai is an even variable and A2, A3 are odd variables, A2 being a column vector 
and A3 a row vector. The fundamental nonvanishing Berezin brackets (here the 
brackets are the ones defined by Eq. (^ in the 2N phase space) are 

{A^{x),p''{y)}^ = 6,''6'{x-y), (55) 

{'^{x),Pi,{y)}j^ = S^{x-y), (56) 

{^{x),p^{y)}^ = 6^{x-y). (57) 

The consistency conditions are 

{(/.I, Hp}j, = dip' - e ^ 7V ^ 0, (58) 

{02, Hp}j^ = -i {d.j - ieAj) ij^^ -eij 7% -m^ +^37° ~ 0, (59) 

{03, Hp}j^ = -I {dj + ieAj) 7^7/; + e7%V^ ^ mi) - 27° A2 ^ 0. (60) 

The last two ones will determine A2 and A3 while the first one will give rise to 
the secondary constraint 

X = dip" — 6-0 7°'0 ~ 0, (61) 

for which the consistency condition will be identically satisfied with the use of 
the expressions for A2 and A3 given by Eq. ( |59D and Eq. (0). Taking the Berezin 
brackets among the constraints we have as nonvanishing results 

)„ , (^3),}, = -i { {i- /)„ , (P,,)^}^ = -» (7%„ . (62) 
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{(02),,x}B = -e{(P^L,^7V}^ = e(^7°)^, (63) 

{(h)^ , x}b = -4{P^)^ ' ^ ^''^}b = -" (^°^). ' ^^^) 

where we explicitly wrote the spinor indexes. Obviously the (j)i constraint is 
first class, but we have another first class constraint. This can be seen from the 
supermatrix A formed by the Berezin brackets among the second class constraints 
X, 02 and 03. Numbering the constraints as $i = x^ ^2 = <p2 and $3 = 03 we 
have this supermatrix in normal form (see ref.|]^, appendix D) given, with spinor 
indexes indicated, by 

{$i,<|.i}^ {^1,(^2) Jb 

A= {($2),,$i}^ {i^2)sA^2)JB 

^{(*3),,$l}^ {(*3),,($2)4^ 

/ 

A = 



$3 







7 



-e [ip 'J 




e (7V) 



($3),} 



(65) 



-z (/; 



/3 
/35 



(66) 



/ 



This supermatrix has one eingevector with null eingevalue that is 

^ 1 ^ 

so there is another first class constraint given by 



V 



(67) 



ip = X + ^e ($2)„ (^j^ - ^e ($3)^ Wp , (68) 

ip = dip' + ie (p^i)+ ^ p^) , (69) 

that we will substitute for x- So, we have the first class constraints 0i and tp, and 
the second class ones 02 and 03. The supermatrix A now reduces to the Berezin 
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brackets among the second class constraints 02 and 03 and is given by 

(70) 

having as inverse 

A-^J " '<^°>-V (71) 

With these result, the Dirac brackets among any variables F and G are 

{F{x),G{y)}^ = {F{x),Giy)}^ 

-^ / dh ({F (x) , ($2), (.)}^ (/)^^ {($3), (z) , G (y)}^ 
+ {F (x) , ($3). {z)}s (/),^ {(<52), W , G (y)})^ . (72) 
The nonvanishing fundamental brackets now will be 

{A,{x),p'{y)}^ = 6,'6'{x-y), (73) 

{Wx i^) , {^)^ (y)}^ = -^ (/),^ S' (x-y), (74) 

{W, (x) , (p^)„ (2/)}^ = 6xJ' (x-y). (75) 

Now we can make the second class constraints as strong equalities and write 
the equations of motion in terms of the Dirac brackets and the extended Hamil- 
tonian given by 

He= f HEd^x = f (nc + Ai0i + a<^) d^x. (76) 

We must remember that, when making 02 = 03 = the constraint ip becomes 
identical to the original secondary constraint x- Then, the equations of motion 
will be 

A ~ {A\ He}^ = -p' + d'Ao - &a, (77) 
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a'^{A',He}^ = Xu (78) 

J9 ~ {p\ He}^ = d,F'' - e ^ tV, (79) 

^K, {^, Ee)j:) = -ieAoip - 7%^' {dj + ieAj) i/; - zm7°'0, (80) 

V^^|z^,ifi5} = ieAo V^ - (Sj - zeAj) ^ 7^7° + im ^ 7°. (81) 

Multiplying Eq. ( PO] ) from the left by 27° we get 

^7° ^= eAo7°?/^ - i7^' {dj + ieAj) ip + mi), (82) 

i {df, + ieA^) -f^ij -mtp = 0, (83) 

while multiplying Eq. (^) from the right by Z7° we get 

« ^ 7O = _e ^ 7OA0 - i (a,- + ieAj) ^p 'y^ - m V', (84) 

« V' T^* (a^. -^eA^) +mij=0. (85) 

These are the equations of motions with full gauge freedom. It can be seen, 
from Eq. (^), that A'^ is an arbitrary (gauge dependent) variable since its time 
derivative is arbitrary. Besides that, Eq. ( [T7| ) shows the gauge dependence of A* 
and, taking the curl of its vector form, leads to the known Maxwell equation 

^ = -i-V(^-«)^^ = -Vxi. (86) 

Writing j^ = e ip 'y^ip we get , from Eq. (^), the inhomogeneous Maxwell 
equation 

^ =V X i - 7, (87) 

while the other inhomogeneous equation 

V ■ E= f (88) 

follows from the secondary constraint (pT]). Expressions (^) and (|53D are the 
known Dirac's equations for the spinor fields ip and ■0- 
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6.2 Hamilton- Jacobi formalism 

Now we apply the formalism presented in the previous sections. From the mo- 
menta definition we have the ^^Hamiltonians" 

Hl, = J{Po + nc)d'x, (89) 

H[ = jp^Sx, (90) 

H!,= f{p^-zij 7°) d'x, (91) 



H', = J p-d'x, (92) 

which are associated, respectively, to t = t*^ (remember that Pq is the momentum 
conjugated to t), Aq, ip and -0. The first two H' are even variables, while the last 
two are odd. Then, using Eq. (pGf ), we have 



dA' = ^dt = f-p' + d'Ao) dt. (93) 

OPi ^ ^ 



From equation (pOD we have 



dp' = -^dt = (djF^' - e i) 7V) dt, (94) 

dpO = --^rft = (a^y -e^j 7V) dt, (95) 

r TTI 

dp^ = — YT^^ ^ i"^^^ ^J-f^ -e^ jf'Af, - m ^) dt, (96) 

dp- = -^^dt + -^^dtP = (-i-i^dj^P + eA^-i^'i) + m^) dt - i7°#. (97) 
The integrability conditions require dH' = 0, which implies for H[ 

dH[ = d fpod^x = 0^dpo = djp> - e ?^ 7V = 0, (98) 

where we made use of Eq. (|95D . This expression is equivalent to the secondary 
constraint (j^ and has to satisfy 

dH'^ = 0; H'^= f (djp> - e i) 7V) d^x, (99) 
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which is indeed identically satisfied. 



For H'2 we have 



dH2 = d f [p^-i^ 7°) d^x = 0^ dp4, = i(d^) 7°, (100) 

which can't be written as an expression like H' = due to the presence of two 
difi^erentials {dp^ and d tj)) but, substituting in Eq. (^), we get 

i (^d i)) 7° = (-idj i>-f^ -e^ YA^ -m^) dt, (101) 

i.e. 

- i ?/i 7° - idj tp-f^ -e'lp jf'Af, -m^=0, (102) 

^ V' 7^ (a^ -^eA^ + m ^= 0. (103) 

For H'^ we have 

dH's = d f p-d^x = 0^dp- = 0, (104) 

that, similarly to the case above, can be used in Eq. (WJ) giving 

(^-ij^djip + eAf,-f>'ifj + m^) dt - i-^^dip = 0, (105) 

i {d^ + ieA^) 7^^ -mij = 0. (106) 

Finally, we can verify, using the above results, that di^g = is identically 
satisfied. 

Equations (|103|) and (|106|) are identical, respectively, to equations (|^) and 



(p3D obtained in Hamiltonian formalism. Besides that, equations (p4\ } and (|79D 
are equivalent, while Eq. (p3| ) corresponds to Eq. (^) except for an arbitrary 
gauge factor. 
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7 Conclusions 

In this work we presented a formal generalization of the Hamilton- Jacobi for- 
malism for singular systems with Berezinian variables, obtaining the equations 
of motion as total differential equations ( p6|) and (|30D . In these equations, the 



coordinates g" = t" (a = 1, ..., i?), whose momenta are constrained, play the role 
of evolution parameters of the system, together with the time parameter t = t°. 
So, the system's evolution is described by contact transformations generated by 
the "Hamiltonians" H'^ and parametrized by t" (with a = 0, 1, ..., R), were H^ is 
related to the canonical Hamiltonian by Eq. ( P^ and the other H'^ {a = 1, ..., R) 
are the constraints given by Eq. (^). This evolution is considered as being always 
restricted to the constraints surface in phase space, there is no complete phase 
space treatment that is latter reduced to the constraints surface, as in Dirac's 
formalism with the use of weak equalities. 

We should observe that, in the case of systems composed exclusively by even 
variables, all parities are equal to zero and equations (pUf), (|30D, (|3l| ) reduce to 
the results obtained in ref.[^. Furthermore, if the system is nonsingular, we have 
H'^ = except for a = 0, so the total differential equations (p6D and (|30|) will be 
reduced to the expressions given by Eq. (^). 

The integrability conditions (which relation to the consistency conditions in 
Dirac's formalism is discussed in the appendix) were shown to be equivalent to the 
necessity of the vanishing of the variation of each H'^ {a = 0, ..., R), i.e. dH'^ = 0. 

The example presented was chosen for its completeness: it is a singular system 
with even and odd variables and its Hamiltonian treatment contains all kinds of 
constraints (primary and secondary, first and second class ones). This example 
is very illustrative, since it allows a comparison between all features of Dirac's 
and Hamilton- Jacobi formalisms. For example, the fact that the integrability 
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conditions dH!^ = and dH'^ = give expressions involving some differentials dt"" 
is related to the fact that the corresponding Hamiltonian constraints 02 and 03 
are second class constraints and determine some of the arbitrary parameters in 
the primary Hamiltonian (^^. Similarly, the fact that the condition dH[ = 
generated an expression like H'^ = is related to the fact that the corresponding 
Hamiltonian constraint 0i is a first class one (see appendix). 

Finally, we must call attention to the presence of arbitrary variables in some 
of the Hamiltonian equations of motion due to the fact that we have gauge de- 
pendent variables and we have not made any gauge fixing. This does not occur 
in Hamilton- Jacobi formalism since it provides a gauge-independent description 
of the systems evolution due to the fact that the Hamilton- Jacobi function S 
contains all the solutions that are related by gauge transformations. 

8 Appendix: Equivalence among consistency 
and integrability conditions 

In this appendix we will show the equivalence among the integrability conditions 
of the formalism showed above and the consistency conditions in Dirac's Hamil- 
tonian formalism, in a similar way to what was made for usual variables |l^. In 



the notation used in this paper the Dirac's primary constraints are written, from 
Eq. (0), as 

H'^ = P^ + H^{q';Pa)^0, (107) 

where a = 1, ..., R; i = 1, ..., N. The canonical Hamiltonian is given by Hq in Eq. 
(pOD , so the primary Hamiltonian Hp is 

Hp = Ho + Kv'', (108) 
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where the v" are unknown coeficients related to the undetermined velocities Q 
[|]. The ordering of the v°' with respect to the H'^ is a matter of choice, since it 
will simply produce a change of sign, but the natural procedure, that identifies 
v°' and Q , suggests the ordering above as a consequence of the ordering adopted 
in the Hamiltonian (|^). This ordering is also the most natural choice to our 
purpose but is, of course, irrelevant for systems containing only usual variables. 
The consistency conditions, which demand that the constraints preserved by time 
evolition, are written as 

ij> {//;, Hp}^ = {i/;, Ho}^ + {h^, h'^}^ ?"= 0, (109) 

where a,fi = 1,...,R and the Berezin brackets here are that given in Eq. (^ 
defined in the usal 2N dimensional phase space and we have made the explicity 
identification q = v°'. 

Multipling the above equation by dt we get 

dH'^^{H'^,Ho}jt+{H'^,H'^}je = 0, (110) 

where, as before, q" = t" but we are still making a = 1,...,R. At this point 
we can already see that, when Dirac's consistency conditions are satisfied we 
have dH' = satisfied. We must see now that we have dH^ = when Dirac's 
consistency conditions are satisfied. The Hamiltonian equation of motion for Hq 
is 

h'o^ {K Hp}^ = {H'„ H',}^ + {Hi H'J^ q\ (111) 

which, multiplied by dt becomes 

dR'^ ^ {Hi Ho}^ dt + {Hi H'J^ dt". (112) 

Remembering that the "momentum" Pq in Hq is independent of the canonical 
variables g* and pi, we have 

dH'o ^ {Ho, Ho}^ dt + {Ho, H'J^ de = {Ho, H'J^ dr. (113) 
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But, if Dirac's consistency conditions are satisfied, we must have only pri- 
mary first class constraints, otherwise we would have conditions imposed on the 
unknown velocities 1 . So, the preservation of constraints in time will reduce to 

{i7;,i7o}^^o, (114) 



and the right side of Eq. (|113|) will be zero. This is simply a consequence of the 



fact that, once all Dirac's conditions are satisfied, the Hamiltonian is preserved. 
So the condition dH^ = is satisfied when Dirac's consistency conditions are 
satisfied. 

This shows that the integrability conditions in Hamilton- Jacobi formulation 
will be satisfied when Dirac's consistency are satisfied. Similarly, we can consider 
that we have the integrability conditions satisfied so that cIHq = dH' = and 
then Eq. ( |110| ), which is equivalent to Eq. ( p.09| ), imphes that Dirac's conditions 
are satisfied. So, both conditions are equivalent. 

Now, we will consider that these conditions are not initialy satisfied. When 
we have only first class constraints in Hamiltonian formalism we will simply get 
a new constraint from some of the conditions ( |1Q9| ). From Eq. (|110|) we see 
that this will imply an expression like dH' = H'dt {H' ?^ is the secondary 
Hamiltonian constraint) which means that there will be a new H' in Hamilton- 
Jacobi formalism that have to satisfy dH' = 0. 

If we have some second class Hamiltonian constraints the consistency condi- 
tions ( |109| ) will imply a condition over some of the velocities Q . From Eq. ( |110| ) 
we see that, in Hamilton- Jacobi approach, there will be conditions imposed on 
some differentials dt"'. 

Such correspondence among the formalisms can be clearly seen in the example 
presented in this paper. 
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Besides that, Eq. ( 110 ) and Eq. (|112|) can be written as 

dH'^^{H'^,H',}^dr, (115) 

were now a, fi = 0,1, ..., R and the Berezin bracket is again defined in the 2N + 2 
dimensional phase space containing t° and Pq. This equation is obviously identical 
to Eq. (1^), that leads to the integrability condition dH' = 0, and its right hand 
side was showed to correspond to Dirac's consistency conditions. Consequently, 
this expression shows directly the relation among consistency and integrability 
conditions. 

It's important to notice that here we are not considering any explicity de- 
pendence on time, neither of the constraints nor of the canonical Hamiltonian, 
because it is an usual procedure in Hamiltonian approach. But the equations 
of Hamilton- Jacobi formalism were obtained without considering this condition 
and, consequently, remain valid if we consider systems with Lagrangians that are 
explicitily time dependent. 

But Hamiltonian approach is also applicable to such systems (see reference [^, 
page 229) and in this case we can follow a procedure similar to that one showed 
here and demonstrate the correspondence among Dirac's consistency conditions 
and integrability conditions. 

Finally, some words about the simpletic structure. Using Eq. (|38|), we can 
writte Eq. ( p6D and Eq. ( |30D in terms of left derivatives as 

dq' = ^^de, (116) 

dp, = -{-if^^^de, (117) 

where, as before, i = 0,1,. ..,N and a = 0,1,. ..,R. These expressions can be 
compactly written as 

dr]' = E'-'^^dr, (118) 
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were we used the notation 

V'' = q\ V''=P^^V' = {q\Pi), / = (C=1,2; ^ = 0,1,...,N), (119) 

E'' = 5^- {6'^,5\ - (-l)^« S\5'^2) , / = (C; ^) , J = (^; j) ; (120) 

that was introduced in page 76 of reference ||^. The Berezin brackets defined in 
Eq. (|) can be written as 

This simpletic notation allows us to obtain the expression for the total differ- 
ential for any function F (t^ , q"" , p"") in a more direct way. Using it in Eq. (^3|) 
we get 

dF = ^drj^, (122) 



where the use of Eq. ( |118[ ) gives 



^^ = 1^^"!^^'" = ^^' ^^^^ ^'"' ^^23) 



in agreement with Eq. 
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